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C3 1 INTRODUCTION 

> - 

, PDS 456 is a luminous, butlow redshift {z = 0.184) radio- 
• quiet quasar identified in 1997 (Torres et al. 1997). The op- 
$_( tical and infra-red spectra (Simpson et al. 1999) show broad 
, Balmer and Paschen lines (e.g. H/3 FWHM 3000 km s~^), 
strong Fe ll, a hard (de-reddened) optical continuum [fi, oc 
j^-o.i±o.i-j^ and one of the strongest 'big blue bumps' of any 
AGN (Simpson et al. 1999, Reeves et al. 2000). It is also 
radio-quiet (Fscijz = 8mJy; Reeves et al. 2000), and is pre- 
sumably not jet dominated or strongly beamed. PDS 456 
has a de-reddened, absolute blue magnitude of Ms ~ —27 
(Simpson et al. 1999), making it as luminous as the radio- 
loud quasar 3C 273 (z = 0.158, Ms « -26). Indeed PDS 
456 is the most luminous known AGN in the local Universe 
(z < 0.3), its luminosity being more typical of quasars at 
z=2-3, at the peak of the quasar luminosity function. 

PDS 456 was first detected as the X-ray source 
RXS J172819.3-141600 in the ROSAT All Sky Survey (Vo- 
ges et al. 1999). Subsequent ASCA and 7?XTi? observations 
of PDS 456 showed that it was highly X-ray variable (see 
Reeves et al. 2000). In particular, during the 7?Jfr£; observa- 
tion, an X-ray flare occurred with a doubling time of just 15 
ksec, implying that the X-ray emitting region was extremely 
compact, less than 2 Schwarzschild radii (or 2Rs) in size. 
Such rapid variability is very unusual for luminous quasars, 



ABSTRACT 

We present evidence from BeppoSAX and XMM-Newton of extreme X-ray vari- 
ability in the high luminosity radio-quiet quasar PDS 456, the most luminous known 
AGN at z < 0.3. Repeated X-ray flaring is found in PDS 456, over the duration of 
the 340 ksec long BeppoSAX observation. The X-ray flux doubles in just 30 ksec, 
whilst the total energy output of the flaring events is as high as 10^^ erg. Under the 
assumption of isotropic emission at the Eddington limit, this implies that the size of 
the X-ray emitting region in PDS 456 is less than 3 Schwarzschild radii, for a lO^M© 
black hole. From the rates of change of luminosity observed during the X-ray flares, we 
calculate lower limits for the radiative efficiency limit between 0.06 and 0.41, implying 
that accretion onto a Kerr black hole is likely in PDS 456. We suggest that the rapid 
variability is from X-ray flares produced through magnetic reconnection above the disc 
and calculate that the energetics and timescale of the flares are plausible if the quasar 
is accreting near to the maximum Eddington rate. A similar mechanism may account 
for the extreme rapid X-ray variability observed in many Narrow Line Seyfert Is. In 
the case of PDS 456, we show that the X-ray flaring could be reproduced through a 
self-induced cascade of ^ 1000 individual flares over a timescale of the order 1 day. 
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as the variability timescale is thought to increase with lu- 
minosity, and black hole mass (e.g. Turner et al. 1999) . One 
possibility is that the accretion rate is unusually high in 
PDS 456, perhaps close to Eddington. An analogy might 
then be drawn with the extreme events observed in several 
Narrow Line Seyfert 1 galaxies (e.g. Boiler, Brandt & Fink 
1996, Leighly et al. 1999), thought to have smaller black hole 
masses (lO^'M© - lO^M©), accreting near to the Eddington 
rate. We report here on X-ray observations of PDS 456, con- 
ducted with BeppoSAX and XMM-Newton in February and 
March 2001. The prime motivation was to study the extra- 
ordinary variability of PDS 456 with an imaging X-ray tele- 
scope, thus negating the possibility of source contamination 
which may occur within the field of view of a non-imaging 
instrument. 



2 THE X-RAY OBSERVATIONS 

PDS 456 was observed by BeppoSAX between 26th February 
2001 and 3rd March 2001, with a total duration of 345 ksec. 
Lightcurves were extracted from circular regions of 4' and 
6' radius around PDS 456, for the MEGS (Medium En- 
ergy Concentrator Spectrometer; Boella et al. 1997) and the 
LEGS (Low Energy Concentrator Spectrometer; Parmar et 
al. 1997) detectors respectively. Background lightcurves were 
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Figure 1. Background subtracted BeppoSAX lightcurves for the 
(a) LEGS (0.3-2 keV) and (b) MEGS (1-10 keV), binned into or- 
bital bins (96 minutes). Repeated X-ray flaring is seen throughout 
the observation. The largest changes are by a factor of xl.9 in the 
MEGS in 35 ksec and a factor of x4.1 in the LEGS in 40 ksec, 
after 230 ksec. The variations imply that for PDS 456, with a 
IO^Mq black hole, the X-ray emitting region is no larger than 3 
Schwarzschild radii (3ils)- 

taken from circular regions, offset from the source. The av- 
erage net source count rates obtained were (1.045 ±0.001) x 
10-1 ct s-i and (7.25 ± 0.12) x lO'^ct s^^ for the MECS 
and LEGS respectively, whilst the background rates were 
less than 10% of the source count rate. XMM-Newton also 
observed PDS 456 on 26th February 2001, with a duration 
of 40 ksec. The start of the XMM-Newton observation was 
coincident with the onset of the BeppoSAX exposure. Here, 
we show only the timing data from the EPIC-pn detector 
(Struder et al. 2001). A full description of the XMM-Newton 
data and spectroscopic observations will be presented in a 
subsequent paper (Reeves et al. 2002, in preparation). Val- 
ues of Ho — 50 kms^^Mpc^i and go = 0.5 have been as- 
sumed throughout and errors are quoted at the 90% confi- 
dence level. 



3 THE X-RAY VARIABILITY OF PDS 456 

The background subtracted BeppoSAX lightcurves are plot- 
ted in Figure 1, grouped into orbital length bins (96 min- 
utes), and extracted over the energy bands 0.3-2 keV and 1- 
10 keV, for the LEGS and MECS respectively. Large changes 
in flux are seen throughout the BeppoSAX observation. The 
most extreme events are the flares seen toward the end of 
the observation at 230 ksec, where the MECS count rate 
increases by a factor of xl.9 over 35 ksec, whilst the LECS 
flux increases by x4.1 over 40 ksec. 

The XMM-Newton EPIC-pn lightcurve is plotted in 
Figure 2 (upper panel), extracted from 0.3-10 keV. The in- 
crease in flux corresponds to the initial rise of the BeppoSAX 
lightcurves over its first 40 ksec of observation, where these 
two observations were concurrent. There is no variability 
over timescales as short as ~ 1 ksec, consistent with a large 
(lO^Mo) black hole mass in this quasar. Using the XMM- 



Newton data, we also searched for evidence of spectral vari- 
ability in PDS 456. A 'softness ratio' was defined, taking 
the ratio of the pn count rate over the 0.3-1.0 keV and 2-10 
keV bands; the time-averaged value was then renormalized 
to 1. A plot of softness ratio against time is shown in Figure 
2 (lower panel). The spectrum becomes softer (by a factor 
~20%) when the count rate is higher during the course of 
the flare in the pn lightcurve. This is also illustrated in Fig- 
ure 3, which plots the softness ratio as a function of total 
pn count rate, the correlation is significant at > 99.99 % 
confidence using a Spearman-Rank test. 

We also constructed cross correlation functions in or- 
der to search for any soft to hard time lags using both the 
BeppoSAX and XMM-Newton data. No delays were found, 
the upper-limit was < 1 ksec. It appears that both the soft 
(0.3-1.0) keV and hard (2-10 keV) bands vary coherently on 
timescales much shorter than that of the overall flare dura- 
tion. For simple reprocessing models, where the thermal disc 
emission is Compton up-scattered to reproduce the hard X- 
ray power-law (e.g. Czerny & Elvis 1987), this implies that 
the size of the reprocessing region is < 10^^ cm or < Q.lRs 
(for a Thomson depth of t ~ 1). The spectral softening 
appears consistent with observations of some AGN, where 
the X-ray spectra are generally softer at higher fluxes (e.g. 
Vaughan & Edelson 2001). 



4 THE SIZE AND EFFICIENCY OF THE 
X-RAY EMISSION IN PDS 456 

Using light-crossing arguments, and assuming that relativis- 
tic beaming is unimportant, one can calculate the overall size 
of the X-ray emitting region in PDS 456 from the expres- 
sion / = ct/{l + r), where t is the rise-time of the flares and 
r is the Thomson depth of the X-ray emitting region. Past 
multi-wavelength studies have shown that PDS 456 has a 
total bolometric luminosity of 10*^ erg , peaking in the 
optical- UV band (Simpson et al. 1999, Reeves et al. 2000). 
With the assumption of isotropic emission (noting that PDS 
456 is a radio-quiet, non-blazar AGN), this then requires 
a 10^ Mg black hole accreting at the Eddington rate. As 
the X-ray flux from PDS 456 can double on timescales of 
~ 30 ksec, then the X-ray emitting region has a maximum 
size oil < 10^® cm, or < 3iis for a IO^Mq black hole, consis- 
tent with the previous RXTE measurement in 1998 (Reeves 
et al. 2000). If one relaxes the assumption that PDS 456 
is accreting at the Eddington rate, and thus the black hole 
mass is > lO^M©, or if the Thomson depth of the region is 
large (r > 1), then the restriction on the radius of the emit- 
ting region is even more severe (i.e. I « 3Rs)- Thus the 
extreme variability suggests that the accretion rate in PDS 
456 is indeed very high, with the X-ray emission originating 
from a very compact region (of < 3Rs), presumably close 
to the event horizon of the putative W^Mq black hole. 

It is possible to calculate the efficiency of converting rest 
mass into energy from the ffaring events. Firstly we estimate 
the total X-ray luminosity from the XMM-Newton spectral 
fits to PDS 456. The XMM-Newton EPIC MOS and pn spec- 
trum of PDS 456 is shown in Figure 4; extrapolating a simple 
power-law flt to the 2-10 keV spectrum shows a clear excess 
of counts below 1 keV in both detectors. Fitting the con- 
tinuum with a hard power-law (of T = 1.98 ± 0.02) and a 
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Figure 2. XMM-Newton EPIC-PN lightcurvc of PDS 456, con- 
current with tlie first 40 ksec of BeppoSAX observation (up- 
per panel). The observations arc consistent with the BeppoSAX 
Ughtcurves, an increase in flux is seen after 20 ksec in XMM- 
Newton but there axe no rapid changes over shorter timescales. 
The lower panel shows the variation in softness ratio (0.3-1.0/2- 
10 keV count rate) versus time, which shows the source softening 
towards the end of the observation. 
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Figure 3. Softness ratio verus total PN count rate. A significant 
correlation is seen, showing that the spectrum of PDS 456 softens 
with increasing flux. 



single temperature blackbody component to model the soft 
excess (with kT = 101 ± 5 eV), then the time-averaged 0.3- 
10 keV luminosity obtained is 1.10 x lO''® erg s~^, compared 
to the bolometic output of 10'*'^ erg s~^ (Simpson ct al. 1999, 
Reeves et al. 2000). Note that slightly higher values for the 
luminosity are obtained if one uses either a Comptonised 
blackbody or a disk blackbody spectrum. 

The time-averaged count rate in the EPIC-pn detec- 
tor is 3.33 ± 0.01 ct s~^, hence one can calculate a con- 
stant factor for the conversion of count rate to luminosity of 
fpn = (3.30±0.01) X 10"'' erg ct~^ From a linear fit to the in- 
crease in the XMM-Newton lightcurve after 20 ksec (Figure 
2), we find a change in PN count rate of 1.64 ± 0.13 ct s~^, 
over 18/(1 + 2:) ksec (quasar rest frame), a factor of xl.62 in- 
crease. Thus the corresponding rate of change in luminosity 
of PDS 456 is AL/At = (3.56 ± 0.27) x 10*^ erg s'^ As- 
suming photon diffusion through a spherical mass of accret- 



ing matter in which the opacity is dominated by Thomson 
scattering, the observed change in luminosity implies a min- 
imum efficiency of ?? > (AL/At)/ {2 x 10*^) (Fabian 1979, 
Guilbert, Fabian & Rees 1983). The derived efficiency, from 
this expression, of > 0.18 ± 0.02, exceeds the theoretical 
maximum for a Schwarzschild black hole, but is consistent 
with the limits for a Kerr metric (Thorne 1974). 

Constraints can also be placed on 77 from the Bep- 
poSAX LEGS lightcurvc. A 0.3-2.0 keV luminosity of 9.65 x 
10"' erg s^^ was derived from a power- law plus black- 
body fit to the (high signal to noise) XMM-Newton spec- 
trum. The mean LEGS count rate, over the portion of 
the lightcurve simultaneous with XMM-Newton was (6.5 ± 
0.2) X lO"'^ ct s~^. Hence the constant factor for convert- 
ing between luminosity and BeppoSAX LEGS count rate 
is fLECS = (1.48 ± 0.07) X lO""" erg ct"^ The two fastest 
events observed here correspond to an increase in count rate 
of (5.42 ± 1.53) X 10-2 ct s'^ in 11500/(1 -I- z) sec after 
230 ksec, and an increase of (4.98 ± 0.95) x 10"^ ct s"^ 
in 17280/(1 + z) sec (after 253 ksec), corresponding to 
AL/At = (8.3 ± 2.3) x 10"*^ erg s'^ and AL/At = (5.1 ± 
0.9) X lO*'^ erg s"^ respectively. The implied efficiency fac- 
tors are then r] > 0.41 ± 0.11 and rj > 0.26 ± 0.05, near 
the maximum permitted value of 77 ~ 0.3 for extraction of 
energy around a Kerr black hole. 

As discussed by Brandt et al. (1999), there are several 
caveats to note about the standard definition of the effi- 
ciency limit (Fabian et al. 1999, Guilbert et al. 1983). The 
limit can be invalid if the X-ray emission is relativistically 
boosted, as can occur in a blazar like AGN, although we 
note here that PDS 456 is radio-quiet (the radio-loudness 
of PDS 456, defined as the ratio of 5 GHz to B band flux, 
is F(icm/FB ~ 0.18) - hence any X-ray emission associated 
with a radio-jet is likely to be weak or absent. Mild relativis- 
tic boosting may be possible through coronal flare emission, 
in this case velocities of /3 = v/c = 0.3 can be achieved 
(e.g. Beloborodov 1999). However in PDS 456 the relatively 
steep (r = 2.0) time-averaged 2-10 keV photon index, and 
the spectral steepening observed during a flare (Figures 2 
and 3), makes this possibility unlikely. 

The derived efficiency also assumes that the radiation 
is emitted at the centre of a spherical region; the limit can 
become invalid if the X-ray emission occurs in the outer 
layers of a sphere of high (t >> 1) optical depth, as the 
flare rise-time can then become unusually short (see the dis- 
cussion in Brandt et al. 1999, Appendix A). Such a high 
value for r seems unlikely in PDS 456, as from light-crossing 
arguements this would require a very compact region of 
I « 3Rs for the X-ray emission. Nonetheless, in order to 
adopt a more conservative approach in estimating rj, we cal- 
culated the total integrated increase in X-ray emission dur- 
ing the course of a flare, rather than using the rise-time of 
the event. The integrated energy emitted during the largest 
flare during the BeppoSAX observation (after 230 ksec) was 
(9.1 ± 0.6) X 10^° erg, over a duration of 103.6/(1 -t- z) ksec. 
Hence the rate of change of X-ray luminosity over the whole 
flare is AL/At = (1.19 ± 0.08) x 10^^ erg s"^ and the de- 
rived efficiency limit is then 77 > 0.060 ± 0.004, close to the 
maximum possible value for accretion onto a Schwarzschild 
black hole. 



© 0000 RAS, MNRAS 000, 000-000 



4 



Reeves et al. 




channel energy (keV) 



Figure 4. The XMM-Newton EPIC MOS (black) and pn (grey) 
spectrum of PDS 456. A power-law of photon index F = 1.98 ± 
0.02 has been fitted to the 2-10 kcV spectrum. Extraploation of 
this power-law model to lower energies (0.3 keV) results in a clear 
soft X-ray excess, which can be modelled by either blackbody 
or comptonised blackbody emission. The total X-ray luminosity 
measured from 0.3-10 keV is then 1.1 x 10*® erg s~^. 



5 DISCUSSION 

The tight limits placed on the timescale and size of the X- 
ray variations in PDS 456 allow us to place constraints on 
its physical origin. The rapid rise times and large ampli- 
tudes of the X-ray variations suggest that those are coherent 
events from the innermost regions {Rin ^ 3Rs) of the accre- 
tion disc. Moreover the high value of the radiative efficiency, 
measured between rj > 0.06 (conservative) and r) > 0.41 (op- 
timal), is indicative of accretion onto a Kerr black hole. The 
timescale of the flares (ifiarc ^ 100 ksec) is shorter than 
the dynamical time scale (tdyn) of the disc at 3Rs (here 
idyn ~ 350 ksec at 3Rs, for a pseudo-Newtonian potential, 
Abramowicz 1996). We suggest that X-ray flares produced 
via magnetic reconnection events in the corona of the in- 
ner disc are a plausible mechanism to explain the observed 
variations over such a short (sub-orbital) time scale. 

Magnetic reconnection events are likely to arise as a 
direct result of buoyant magnetic flux tubes emerging ran- 
domly from the disc surface. Reconnection then takes place 
in regions of the corona in which oppositely directed mag- 
netic field lines come in close contact. Di Matteo (1998) ex- 
amined coronal heating and flare production via Petschek 
reconnection events in AGN. The discussion presented here 
is based on many of the results of that analysis. Flare events 
in many AGN can be accounted for by the onset of an 
ion-acoustic instability associated with slow MHD shocks 
and Petschek reconnection in the accretion disc corona. The 
fundamental energy source driving this process is the or- 
bital energy of the inner disc, as shear stresses in the turbu- 
lent disc can rapidly amplify any seed fields. The resultant 
flux tubes are subject to a buoyancy (Parker) instability 
and emerge from the disc into the magnetically-structured 
corona. In this way a fraction of the accretion energy is con- 
verted into coronal magnetic energy. The emergence rate 
of magnetic energy into the disc corona should therefore be 
proportional to the the kinetic energy flux through the inner 



disc -Edisc ~ dt{2-KRT,v'^) ~ Mt;^ ~ M(? jt in steady state 
at i? = ei?s, where E is the disc surface density, is the 
azimuthal velocity and M is the steady state accretion rate. 
Assuming a 10^ M© black hole accreting at the Eddington 
rate, we have _Edisc(3-Rs) ~ lO*'^ erg s~^ with e = 3 and from 
observation the total output of the flares is iSflare lO^'^ erg, 
over the lifetime (recurrence time) of the flares of tree ~ 1 
day. A relatively high efficiency (> 0.1) for the conversion 
of accretion energy to X-ray luminosity is required to power 
the flares in PDS 456, and is another indication that the 
massive black hole in this object is accreting at a rate close 
to the Eddington limit. 

The buoyancy condition for magnetic flux tubes within 
the accretion disc is that the magnetic pressure should ex- 
ceed the local gas pressure (e.g. Coroniti 1981), the critical 

condition being ~ ridiscfeTdisc, where Be is the local mag- 
netic field strength, ridisc is the particle density within the 
disc and Tdisc is the disc temperature. The structure of the 
accretion fiow close to the last stable orbit is unknown, how- 
ever estimates can be made by assuming that accretion takes 
place via a thin disc. Using a steady-state Shakura-Sunyaev 
disc with a viscosity parameter a ---^ 0.1 yields the estimates 
Tidisc ~ 10^"^ cnr-^ and Tdi^c - lO'^ K at r ~ 3i?s. Applying 
these estimates to the magnetic pressure condition above, 
leads to a value for the field strength of a buoyant flux tube 
within the inner regions of the disc of Be ^ 5000 G. 

The observations provide two constraints on the den- 
sity of the X-ray emitting region. These arise from the con- 
straints on the fitted column density and the ionisation pa- 
rameter of the highly ionised absorber in PDS 456. This is 
apparent in the hard X-ray spectrum of PDS 456 above 7 
keV, in the form of deep K-shell edges of highly ionised iron 
(Fe XXV and Fe xxvi), present both in the current XMM- 
Newton EPIC spectrum (Reeves et al. 2002, in preparation) 
and in earlier ASCA and RXTE observations (Reeves et al. 
2000). From the XMM-Newton observations, the X-ray spec- 
trum of PDS 456 is well fit with a column density nn ~ nRin 
in the range 9 x 10^^ - 5 x lO^'' cm"^ implying 9 x 10* cm"* 
< n < 5 X 10® cm^"'. On the other hand, the ionisation pa- 
rameter U ~ L^/nR^-^ lies in the range 3 x 10** - 2.5 x 10*, 
where Lx ~ 10** erg s^^ is the hard X-ray ionising fiux 
leading to the limits 4 x 10^ < n < 3 x lO" cm"*. If we 
associate this density with the hot post shock gas we obtain 
the estimate ~ 5 x 10^ cm~*. The evolution of the buoy- 
ant fiux tubes provides a natural explanation for tins very 
low value of the density of the X-ray emitting region. As the 
flux tubes emerge from the disc and rise into the corona, 
gravitational downflow of plasma causes the density within 
the flux tube to decrease to the point at which reconnection 
can take place efficiently. The density estimate above implies 
that the buoyant flux tubes rise far into the corona before 
reconnection occurs. An estimate of the height at which re- 
connection occurs in terms of the pressure scale height of 
the disc {H) can be obtained from (Di Matteo 1998) 

/ N 1/4 

/Iflare-^f^j H ^ ^QH . (1) 

In the case of a Shakura-Sunyaev accretion disc H/R < 
0.1 implying that /iflarc ^ Rin- Hence the total magnetic 
energy stored in the corona of the inner disc should be of 
the order B^RfjSn ~ 10*^ erg. This large reservoir of stored 
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accretion energy is enough to explain the large amplitude X- 
ray variations seen in PDS 456, and highlights the efficiency 
with which accretion energy must be converted to magnetic 
energy to explain these flaring events. 

Wc now consider the rise time of the events in the 
light curves (trise ^ 30 ksec) in the light of the above 
model. An attractive feature of the Petschek model is the 
short time scale on which reconnection can take place. This 
can be of the order of a few Alfven times ta ~ 1/va ~ 
Rin{4:'!vnsmp)^^^ /Be whore / is the characteristic length scale 
associated with changes of the magnetic field (here I ~ Rin), 
va is the local Alfven speed and rup is the proton mass. The 
estimates above yield ta ~ 50 ksec, encouragingly close to 
the observed rise time of the flares. The time scale ta pro- 
vides a constraint on the density of the accretion disc and 
post shock gas. An accretion disc density significantly lower 
than ndisc ~ lO^'^ cm~^ or a shock density much higher than 
Us ~ 10^° cm^^ would increase ta to a level inconsistent 
with the observations. 

The energy production rate for an individual flare can 
be estimated from (Di Matteo 1998):- 

-C'flare ~ j/j ~ 1/2 l^J 

me me 
where L is the length of the slow shock region, Tflare is the X- 
ray temperature and rUe is the electron mass. The dimension 
L is the essentially the length of region of oppositely-directed 
magnetic field lines which wc constrain to be of the order 

L ~ i?in. Adopting Tflare ^ 10^ K we obtain i;flare < 10"*^ 

erg s~^ . This suggests that the large scale X-ray variations in 
the light curve of PDS 456 involve > 1000 or so individual 
flares. These flares could not produce the large amplitude 
variations observed in PDS 456 if they were incoherent. 

However it is possible that the magnetic structure 
within the disc corona could reach a self-organized, criti- 
cal state in which the reconnection and flaring of one flux 
tube could prompt similar flares in its neighbors, allowing 
a coherent cascade of flares to develop (Leighly & O'Brien 
1997). The timescale on which the global cascade could take 
place can be constrained from the Alfven time in the region 
involved, which is consistent with observed timescale of the 
flaring events in the light curve. 

This leads to a suggestion of why it is that PDS 456 
releases such a large fraction of its accretion energy in the 
form of large-amplitude, coherent X-ray variations. A flare 
cascade along the lines of that suggested above would re- 
quire that magnetic energy was stored in the disc corona 
until some triggering criterion was reached, promoting the 
flrst flare event. The energy involved in such a cascade is 
likely to be governed by the time scale on which magnetic 
energy is pumped into the disc corona (tmag) and the time 
scale on which the triggering criterion is satisfied (ttrigger). 
Two possible cases emerge! (a) imag ^trigger in which case 
the flares will be incoherent and the resultant X-ray vari- 
ability small-amplitude and (b) tmag ^ ^trigger producing 
the possibility of large-amplitude coherent flare cascades. 
Unfortunately little is known about ftrigger, however tmag 
must be related to -Bdisc via tmag ~ Ec/E^mc where Ec is 
a measure of the minimum magnetic energy required to be 
stored in the disc corona before it can reach a critical state 
allowing a flare cascade. Hence fmag oc Ec/M oc Ec/mMBH, 
where rh is the accretion rate in Eddington units. Hence 



systems in which a high mass black hole accretes at the Ed- 
dington limit should be more able to release stored magnetic 
energy in the form of a coherent flare cascade, providing at 
least a tentative explanation of why the X-ray variability of 
PDS 456 is so extreme. The similarity of the X-ray behavior 
of PDS 456 and the Narrow Line Seyfert Is (see Leighly et 
al. 1999) suggests that m is the critical factor in determining 
the scale of the X-ray variability in AGN. The energy of X- 
ray flares in some NLSls can be as high as ~ 10** erg (e.g. 
IRAS 13324-3809, Boiler et al. 1997), a factor of - 10'^ loss 
than the flares observed in the light curve of PDS 456, in 
good agreement with the ratio of black hole masses in these 
systems. If we use these observations to make the connection 
Ec oc Mbh, then tmag ~ l/rh. Such a relation would suggest 
that any system accreting close to the Eddington rate would 
be likely to show strong X-ray variability. 



6 CONCLUSIONS 

Recent BeppoSAX and XMM-Newton observations have 
shown that the luminous quasar PDS 456 exhibits rapid X- 
ray variability on timescales of ~ 30 ksec, with a total energy 
output of 10^^ erg s~^ for the flaring events. This limits the 
size of the X-ray emitting region to < 3 Schwarzschild radii 
for a 10^ M0 black hole. The energetics of the of the accre- 
tion disc in PDS 456 can power its extreme X-ray variability 
if the black hole is massive (> 10^ Mq) and is accreting 
close to the Eddington rate. Coronal magnetic flare events 
can explain the X-ray variability as long as the disc is able 
convert accretion energy into coronal magnetic energy effi- 
ciently, and that this energy can be released in the form of 
a self-induced cascade of ^ 1000 individual flare events on 
a timescale of the order 1 day. 
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